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ABSTRACT 

We have extended our previous analysis of morphologically selected elliptical and SO 
galaxies in the Hubble Deep Field (HDF) North to include HST data in the HDF 
South and the HDFS-NICMOS areas. Our final sample amounts to 69 E/SO galax- 
ies with K < 20.15 over an area of 11 square arcmins. Although a moderately small 
number over a modest sky area, this sample benefits of the best imaging and photo- 
metric data available on high-redshift galaxies. Multi- waveband photometry allows us 
to estimate with good accuracy the redshifts for the majority of these galaxies which 
lack a spectroscopic measure. We confirm our previous findings that massive E/SOs 
tend to disappear from flux- limited samples at z > 1.4. This adds to the evidence that 
the rest-frame colours and SEDs of the numerous objects found at 0.8 < z < 1.2 are 
inconsistent with a very high redshift of formation for the bulk of stars, while they are 
better consistent with protracted (either continuous or episodic) star-formation down 
to z < 1 . These results based on high-quality imaging on a small field can be comple- 
mented with data from colour-selected Extremely Red Objects (EROs) on much larger 
sky areas: our claimed demise of E/SOs going from z = 1 to z = 1.5 is paralleled by a 
similarly fast decrease in the areal density of EROs when the colour limit is changed 
from {R — K) = 5 to (R — K) = 6 (corresponding to z ~ 1 and z ~ 1.3 respectively). 
Altogether, the redshift interval from 1 to 2 seems to correspond to a very active phase 
for the assembly of massive E/SO galaxies in the field, and also probably one where a 
substantial fraction of their stars are formed. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: 
formation - galaxies: photometry - infrared: galaxies. 



1 INTRODUCTION 

The evolution with redshift of the mass function of bright 
galaxies is an important discriminant among competing cos- 
mogonic scenarios. Present instrumentation, including HST 
and large optical telescopes on ground, is close to provide 
decisive tests of these alternative schemes for galaxy forma- 
tion. 

The surface density of Extremely Red Objects (EROs) 
has been recently used by Daddi et al. (2000) to constrain 
the epoch of formation and the evolution of elliptical galax- 
ies in the field. These authors claim that the density of 
EROs with (R — K) > 5 (the typical colour of an old el- 
liptical galaxy at z > 1) is consistent with the passive evo- 
lution of giant ellipticals, implying an early formation of 
the spheroidal population through violent starbursts at high 
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redshifts (z > 2.5). This formation paradigm for ellipticals 
has the nice feature to make clear and testable predictions. 

On the other side, hierarchical galaxy formation models 
(Kauffmann & Chariot 1998) predict that elliptical galaxies 
form through the merging of disk galaxies at moderate red- 
shifts. This scenario is supported by other studies empha- 
sizing a lack of massive ellipticals at redshift greater than 
z ~ 1.5 (Zepf 1997; Franceschini et al. 1998, hereafter FA98; 
Barger et al. 1999). 

In particular, FA98 studied the photometric and statis- 
tical properties of a complete sample of early-type galaxies 
with K < 20.15 in the HDF North. The FA98 sample was 
selected on a deep A'-band image from Kitt Peak, by ap- 
plying a rigorous morphological classification scheme based 
on the observed surface brightness profiles. FA98 found that 
the vast majority of bright early-type galaxies in the HDF- 
N are located at z < 1.3 and display colours indicative of a 
fairly wide range of ages. There is a remarkable absence in 
the FA98 sample of objects with z > 1.3. According to the 
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passive evolution scenario, these objects should be easily de- 
tected, due to the luminous star formation phase expected 
to happen at these redshifts. 

To enforce the statistical significance of these results, an 
extension of the analysis to different sky regions would be 
critical. For this reason, we extend in this paper the study to 
the HDF-South WFPC2 and NICMOS fields, by selecting 
a new complete sample of early-type galaxies and following 
the same prescriptions outlined in FA98 to process the data. 

For consistency with previous analyses, we assume 
Ho = 50 Km s _1 Mpc^ 1 and qo = 0.15, with zero cos- 
mological constant A. 



2 SAMPLE SELECTION 

The present analysis makes use of a collection of high quality, 
deep images, well suited for the study of field galaxies and 
covering three different areas: (a) the HDF-North, (b) the 
HDF-S WFPC2 area and the (c) HDF-S NICMOS field. Our 
galaxy catalog includes early-type objects selected in the 
near infrared in such fields. 

2.1 Optical and near infrared imagery 

We defer the reader to FA98 for details on the HDF-N ob- 
servations. The WFPC2 strategy adopted for the HDF-S 
(Williams et al. 1998) was similar to that adopted for the 
HDF-N. In particular, the images were taken in the same 
four bands (F300W, F450W, F606W, F814W), with simi- 
lar total exposure times. For each filter, several exposures, 
obtained in dithering mode, have been combined in a single 
image with pixel scale of ~ 0.04". The sky area covered by 
both HDFN and HDFS is 10.6 arcmin 2 . 

Da Costa et al. (1998) observed the HDF-S WFPC2 
area from optical to near-IR at the ESO 3.5 New Technology 
Telescope (NTT). The optical observations were carried out 
using the SUSI2 camera in binned mode, yielding a scale 
of 0.16 arcsec/pixel. Infrared observations were obtained in 
the J, H and K bands, using the SOFI camera (scale of 0.29 
arcsec/pixel). 

The HDF-S NICMOS images were acquired using cam- 
era 3 with a pixel size of 0".2 and the filters F110W, F160W 
and F222W. The dithered images were processed, drizzled 
and combined giving a final pixel size of 0".075 and covering 
an angular area of ~1 arcmin 2 . 

The HDF-S NICMOS field has also been observed as 
part of the VLT-UT1 Science Verification program, using 
the Test Camera (VLTTC) at the Cassegrain focus and the 
U, B, V, R, I filters. The total field sampled covers 1.5x1.5 
arcmin 2 and the scale of processed images is of 0.092 arc- 
sec/pixel ( Font ana et al. 1999). 

2.2 The photometric filter 

The criteria adopted to ensure a careful selection of the sam- 
ple in the HDF-N, as well as the accurate treatment of pho- 
tometric data in seven passbands (U, B, V, I: from HST 
imaging, by Williams et al. 1996; J, H, K: from KPNO 
imaging, Connolly et al. 1997) have been discussed in FA98 
and Rodighiero et al. (2000). For the other fields we have 
applied the same selection scheme used by FA98. We refer 



the reader to those papers for a thorough description of the 
procedure, while for convenience we summarize the main 
steps below. 

HDF-South- WFPC2 field 
Our primary selection is in the 7f-band, in order to min- 
imize the biases due to the effects of K and evolutionary 
corrections. The sample of galaxies has been extracted from 
the SOFI-A" image, running SExtractor (Bertin & Arnouts 
1996). As in FA98, to determine the limit of completeness in 
the K band (Kl) for inclusion in our sample, we have used 
large simulated fields, including several HDF-like toy galax- 
ies and mimicking the image quality of the if-band image. 
Running SExtractor on these synthetic frames allowed us to 
check its performance, determining the biases and the stan- 
dard deviations of the magnitude estimates as a function 
of the surface brightness. We then derived the corrections 
that must be applied to the SExtractor photometry and the 
limit for completeness Kl = 20.2 (in the standard system). 
From the simulations we derived also an evaluation of the 
K-band limiting surface brightness observable in the field, 
fix ~ 23mag/arcsec 2 . An analogous procedure was used to 
derive, for each object in the selected sample, the corrected 
magnitudes in the J and H bands. We did not correct optical 
HST magnitudes. 

HDF-South-NICMOS field 
The same procedure has been adopted for the HDF-S NIC- 
MOS field. In this case, however, we have chosen the F160W- 
HST band (~ H ) for source selection, because of the higher 
noise in the F222W image (K band). The importance of this 
field concerns its depth in the near-IR (observed from the 
space), that allows to reach deeper limiting magnitudes. In 
particular, from the simulations we derived a magnitude of 
22.05 for the limit of completeness in the H band. As for the 
HDF-S- WFPC2 field, we computed the corrections for the J 
and K bands, and also for the ground based observations in 
the U, B, V, R and I bands. Note that the selection in the 
F160W band does not introduce appreciable bias with re- 
spect to the nominal K band selection adopted in the other 
fields, because the typical spectra for z ~ 1 E/S0 galaxies 
are flat between 1.5 and 2.2 micron and are contributed by 
similarly old stellar populations. 

2.3 The morphological filter 

We have used the Galaxy Authomatic Surface PHOTome- 
try tool (GASPHOT; Pignatelli and Fasano 1999, 2000) to 
perform the morphological analysis of all galaxies brighter 
than the proper limiting IR-magnitudes in the WFPC2 and 
NICMOS fields of the HDF-S. 

This tool operates on CCD images and consists of three 
different steps: i) first, it exploits the basic extraction capa- 
bilities of SExtractor to perform a simultaneous slicing of 
all selected objects with a given step of surface brightness 
(we used m .l), providing luminosity, ellipticity and posi- 
tion angle profiles for each object as a function of the equiv- 
alent radius of the isophote; ii) secondly, it gives an accu- 
rate, smooth representation of the PSF by averaging (with 
a fe-sigma clipping procedure) the luminosity profiles of the 
stars identified in the frame; iii) finally, it uses the general- 
ized de Vaucouleurs law fi r —^o oc (r/r e ) 1/n (where /x r is the 
surface brightness at radius r and r e is the half light radius; 
Sersic 1968, see also Ciotti 1991), convolved with the proper 
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Figure 1. Histograms of the Sersic index n obtained by running 
GASPHOT on a crowded artificial field containing 200 toy galax- 
ies with different magnitudes and radii. Galaxies with exponential 
profiles (fOO objects) are binned with the solid line histogram, 
whereas the dotted line histogram refers to galaxies with r 1 / 4 
profiles. 



Figure 2. Correlation between the Sersic index (n) and the 
de Vaucouleurs morphological type (T) for a sample of galaxies in 
the cluster Abcll 2670. The morphological types were derived by 
visual inspection of the individual galaxies, complemented by lu- 
minosity and geometrical profiles obtained from detailed surface 
photometry (Fasano et al. 2001). 



PSF, to fit the equivalent luminosity profiles of galaxies, 
providing unbiased estimates of the profile parameters. The 
performances of GASPHOT are presently undergoing a sys- 
tematic investigation (Pignatelli and Fasano 2001) by means 
of both extensive simulations of toy galaxies of different mor- 
phological types and comparison with galaxy samples whose 
morphology has been previously estimated by visual inspec- 
tion of images and profiles. Figures 1 and 2 illustrate some 
preliminary results of this analysis. In particular, in Figure 1 
we report the histograms of the Sersic index n obtained by 
running GASPHOT on a crowded artificial field containing 
200 toy galaxies (with different magnitudes and radii) hav- 
ing exponential (100 objects) or r 1 ^ 4 luminosity profiles. In 
Figure 2 we show the correlation between the Sersic index 
(n) and the de Vaucouleurs morphological type (T) for a 
sample of galaxies in the cluster Abell 2670. In this case the 
morphological types were derived by visual inspection of the 
individual galaxies, complemented by luminosity and geo- 
metrical profiles obtained from detailed surface photometry 
(Fasano et al. 2001). From both the above figures the value 
of n given by GASPHOT appears to be a good indicator of 
the 'true! morphological type. 

According to Figure 1, after running GASPHOT on our 
high resolution images, we decided to exclude from the final 
samples the galaxies with n <2 (considered late-type). In 
addition, we rejected the galaxies whose early-type classifi- 



cation was judged dubious at a visual inspection (3 objects 
in the WFPC2 field and 1 object in the NICMOS field). The 
final sample includes the 35 E/S0 galaxies extracted from 
the HDF-N in FA98, 28 galaxies from the WFPC2 field and 
11 galaxies from the NICMOS field, for a total of 74 ob- 
jects. Note that only 6 galaxies in the NICMOS field are 
brighter than the completeness limit for the WFPC2 fields 
(K L = 20.15). 

In Table 1 we report the data of our sample: coordinates 
(at J2000), photometric redshifts and the observed K mag- 
nitudes (the data for HDFN galaxies have been published in 
FA98). 



3 ANALYSIS OF THE BROAD BAND 
SPECTRA 

3.1 Tools for spectral synthesis 

The observed UV-optical-near/IR spectra of the sample of 
E/S0 galaxies presented in this work have been fitted with 
spectrophotometric models describing their integrated emis- 
sion. In order to check the consistency of our approach, we 
decided to compare the results from two different codes: 
GRASIL (Silva et al 1998) and PEGASE (Fioc & Rocca- 
Volmerange 1997). As in FA98, we have considered two mod- 
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els describing the spectral behaviour of elliptical galaxies. 
We adopt a Salpeter IMF with a lower limit Mi = O.15M 
and a Schmidt-type law for the star formation rate (SFR): 
ty(t) — vM g {t), where v is the efficiency and M g (t) is the 
residual mass of gas. A further relevant parameter is the 
timescale t in f a u for the infall of primordial gas. The evo- 
lution patterns for the two models considered here are ob- 
tained with the following choices of the parameters: Model 
1: ti n f a ii = 0.1 Gyr, v = 2 Gyr~ 2 ; Model 2: t infaU = 1 Gyr, 
v = 1.3 Gyr~ 2 . The corresponding SF laws have a maxi- 
mum at galactic ages of 0.3 (Mod 1) and 1.4 Gyr (Mod 2), 
and are truncated at 1 (Mod 1) and 3 Gyr (Mod 2) to mimic 
the onset of galactic winds. 

A x 2 fitting procedure has been adopted for the deter- 
mination of the best solution for each galaxy spectrum. We 
use GRASIL as our reference model and use PEGASE only 
for comparison (the two give equally good fits) . 



3.2 Photometric redshifts and completeness 

An application of spectrophotometric models is the determi- 
nation of photometric redshifts. Our sample in the HDF-S 
lacks a complete spectroscopic coverage (40% of sample ob- 
jects in the HDF-N have measured redshift). We deemed 
useful to compare the GRASIL and PEGASE redshift es- 
timates. There is a slight tendency for GRASIL to predict 
larger values of z, but in general the agreement is quite good 
(~ 10%). 

The distribution in redshift of a flux limited population 
of galaxies provides a constraint on its evolutionary history 
and on the epoch of formation. However, selection effects 
have to be carefully taken into account, in particular the 
surface brightness limit, the cosmological dimming and the 
K-correction. Our selection in the Tf-band minimizes K- and 
evolutionary corrections as a function of redshift. Our eval- 
uation of the effective radii allows us to control the effects 
of the limiting surface brightness (fix ~ 2Smag/arcsec 2 ) 
observable in the if-band image. We found that the cutoff 
in surface brightness in the K image does not affect the se- 
lection above our limit of Kl = 20.15, for redshifts up to 
z ~ 2.5. 



4 RESULTS AND DISCUSSION 
4.1 Redshift distributions 

Figure 3 reports the observed redshift distribution of ellip- 
tical galaxies in the HDF-S (dashed histogram), compared 
with that in the HDF-N (continuous line). The two distribu- 
tions clearly differ: the HDF-N distribution shows a peak at 
z ~ 1, while in the HDF-S the peak is shifted toward lower 
redshifts (z ~ 0.5). These differences can be explained in 
terms of clustering effects in the limited volumes sampled. 
It is worth noticing the global lack of bright spheroids in 
both fields at high redshifts. In particular, early-type galax- 
ies seem to disappear at z > 1.5. Since this conclusion is now 
based on a significantly larger area than used by FA98, it is 
becoming more unlikely that this is due to galaxy clustering 
effects [see further arguments about this in Sect. 4.3]. 

In Figure 4 the observed distribution (continuous line) 
relative to the total area of ~ 11 square arcminutes (HDFN 
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Figure 3. Observed redshift distribution in the HDFN (solid line) 
compared with that in the HDFS (dashed line). 
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Figure 4. Observed total redshift distribution (solid lines), in- 
cluding our sample objects from the HDFN, the HDFS and the 
NICMOS field at the limit magnitude of K L = 20.15. The data 
are compared with predictions based on Model 1 for different red- 
shift of formation zp (dashed lines); from top left to down right: 
z F = 1.5, 2, 3, 4. 
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+ HDFS-WFPC2 + NICMOS), for K < 20.15, is compared 
with the corresponding distributions predicted in the mono- 
lithic collapse scenario (dashed lines) . The four panels report 
the predictions of Model 1 for increasing redshift of forma- 
tion zf- Low values of zf appear to be consistent with the 
observed distribution (zf ~ 1.5 ~ 2), while an early epoch 
of formation for field E/SO seems to be ruled out (see also 
Nakata et al. 1999). In particular, 39 E/SO galaxies would be 
expected at z > 1.5 for zf = 3 and 50 for zf = 4, compared 
with the 2 observed. Our data are more consistent with the 
hierarchical scenario, where large galaxies are predicted to 
form from the merging of smaller units. This finding is in 
agreement with Menanteau et al. (2000), who estimate that 
at z ~ 1 about half the field ellipticals must be undergoing 
episodes of star- formation. 

Among the solutions proposed to explain the demise of 
high-z ellipticals, FA98 consider the bias in the morphologi- 
cal filter against disturbed E/SO's (consequence of merging), 
which may become significant at z > 1. This point has been 
addressed by Rodighiero et al. (2001), who have analyzed, 
in a way similar to that presented in this paper, a comple- 
mentary sample of 52 late-type spiral and irregular galaxies 
selected in the same HDF-N area with total if-magnitudes 
brighter than K=20.47. One of their main conclusions was 
about the fact that only 2 out of 52 galaxies in the sample 
have z > 1.4. The lack of massive galaxies of any morpho- 
logical kinds at high z supports that the morphological se- 
lection of the spheroidal subsample is not responsible for its 
cutoff at z > 1.4. 



4.2 Ages and baryonic masses for the sample 
galaxies 

The observed galaxy SEDs present two different spectral 
behaviours: red objects with spectra typical of old ellipti- 
cals, and blue objects that are flatter at all wavelengths, 
indicating the presence of residual on-going star formation 
and a spread in the ages of stellar populations. Our \ 2 fit- 
ting procedure lead us to the same conclusions as found 
by FA98: model 2 with protracted star-formation seems to 
better reproduce all kinds of spectral behaviour of galaxies 
in the HDF-S, providing a formally acceptable total x 2 of 
~ 290 (7 band data for 28 HDF-S galaxies), corresponding 
to a reduced chi-square of xt ~ 1-5. Model 1 gives accept- 
able fits only for galaxies with red spectra, while the total 
X 2 ~ 400 is unacceptably high (xt ~ 2). This result sup- 
ports the hypothesis that, on average, star formation activ- 
ity in these galaxies occurs over protracted periods (3 Gyr 
in the model), rather than in short-lived single starburst (at 
least when considering our limited set of models). 

We compare in Figure 5 the rest-frame (V—K) and (B— 
J) colours (computed from best-fitting solutions) of early- 
type galaxies as a function of redshift with the predicted 
colours of single stellar populations with solar metallicity. 
The colours for HDF-N ellipticals (open circles) are consis- 
tent with those found for HDF-S objects (filled squares), 
whereas the large spreads observed in both colour distribu- 
tions correspond to a wide range of ages for the galaxies in 
the sample, typically from 1 to 5 Gyr. The figure shows a 
tendency for high redshift ellipticals to be intrinsically bluer, 
in particular in the HDF-S, while the clumped distribution 
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Figure 5. Rest frame (V — K) (upper panel) and (B — J) (lower 
panel) colours of early type field galaxies, compared with pre- 
dicted values for single stellar populations with solar metallicity. 
Filled squares refer to HDFS objects, open circles to the HDFN. 



observed in the HDF-N at z ~ 1 and V - K ~ 2.5 likely 
indicates the presence of clustering at this redshift. 

We plot in Figure 6 the baryonic masses of our sam- 
ple galaxies as a function of redshift, and compare them 
with predictions of model 1 (dotted lines) and model 2 
(solid lines) corresponding to the lowest galaxy mass still 
detectable within the flux limit of the survey (Kl = 20.15) 
as a function of redshift. Different curves in the figure refer 
to various redshifts of formation zp (numbers attached to 
each curve). 

In agreement with the results of Figure 4, only low val- 
ues of zp (< 2) are consistent with the observed mass distri- 
bution. In particular, the predictions of model 1 with zf — 5 
are inconsistent with the observed distribution, half of the 
sample galaxies having estimated masses lower than the pre- 
diction. This model fails to reproduce the broad-band spec- 
tra of each single galaxy. Model 2 with zf = 2 seems to bet- 
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Figure 6. The observed baryonic mass of the sample galaxies at 
Kl = 20.15 is plotted against the redshift. The curves represent 
the predictions of model 1 (dotted lines) and model 2 (solid lines) 
for the lowest galaxy mass still detectable within the flux limit 
of the survey (K^ = 20.15) as a function of redshift for various 
redshift of formation (number attached to each curve) . The error 
bar shows the mean uncertainty on the determination of the mass 
from spectral fitting, when considering the two different models. 



ter reproduce the observations, and provides an appropriate 
match with the lower envelope of the mass distribution. Sev- 
eral galaxies at z < 1 have masses requiring a lower redshift 
of formation (zp = 1). 

We stress that Figure 6 does not imply that ellipticals at 
z > 1.5 do not exist. Rather, according to the prescriptions 
of hierarchical models (Kauffmann & Chariot 1998), they are 
likely to be missed in our flux limited sample due to their 
lower mass (and luminosity). This is strongly suggested by 
the location in Figure 6 of the five objects selected in the 
NICMOS field and fainter than K = 20.15 (three-pointed 
stars). 

Benitez et al. (1999) have detected and identified in 
the NICMOS field 4 early-type galaxies at z > 1 with 
H(AB) < 25 (corresponding to H ST ~ 23.6). Of these, ET2 
(in the Benitez et al. list) is fainter than our limit of com- 
pleteness, ET1 and ET3 are included in our sample, while 
ET4 has been rejected as it is undetected in the VLT images 
and SEDs based on the NICMOS near-IR data cannot be 
satisfactory fitted with any of our galaxy models. 

Benitez et al. suppose that the presence of luminous 
red galaxies at high z in a small field indicates that evolved 
(zf > 5), massive galaxies are common at z > 1. They 
suggest that the absence of distant ellipticals claimed by 
others (FA98, Menanteau et al. 1999, Barger et al. 1999) may 
be due to the depth of the NICMOS data compared with 
ground-based IR data (about 3 magnitudes) or to clustering 
effects. 

The present study is not inconsistent with the results 



reported by Benitez et al.. First, we confirm the presence 
of a substantial number (21 out of 74) of E/SO galaxies at 
z between 1 and 1.4. We also prove the existence of E/SO 
galaxies above z = 1.5. In particular if we include the fainter 
NICMOS data, we find 4 early-types in NICMOS at 1.5 < 
z < 3 with H < 22.05. Our claim is rather that the baryonic 
masses found for these high-z objects (typically few 10 10 Mq 
units) are typically lower than expected on consideration of 
the very massive (M > 10 11 A/0) galaxies detected around 
z = 1 in our areas. The areal density of massive E/SO as a 
function of redshift shows a drastic decrease, not a complete 
disappearance, above z ~ 1.4. 

All this seems to favour an evolutionary scenario where 
the bulk of luminous spheroids are assembled at relative 
recent cosmic epochs (zf < 2), although not as recent as 
supported by Kauffmann & Chariot (1998). 

Obviously, the assembly of stellar populations into the 
dynamically relaxed shapes of spheroids could have occurred 
later than the birth epochs of stars. The distribution of red- 
shift of formation for stellar populations in E/SO has been 
quantified by FA98 in their HDF-N analysis. This showed 
that the bulk of stars formed at cosmic epochs between z = 1 
and z = 3. 

The results of a similar analysis of our combined sam- 
ple is reported in Figure 7. The histogram in the Figure 
shows the mass in stars formed per unit comoving volume 
in the redshift bins. This estimate is based on our best-fit 
photometric model solutions for all galaxies in the sample, 
by adding up the individual rates of star formation in the 
redshift bins. The mass is corrected for the portion of the lu- 
minosity function not sampled by our K-Qxa. selection using 
the local observed LF in the Ji'-band (see FA98 for details 
on the computation). 

As shown by the figure, the history of formation of stars 
for this combined HDFN, HDFS, NICMOS sample is differ- 
ent from the partial result based on the HDF-N. The dashed 
curve is the normalized integral distribution: in the present 
case and for our adopted cosmology , 30% of stars are formed 
at at z < 1, 40% between z — 1 and 2, and the remaining 
30% at 2 > 2. 

4.3 Effects of cosmic variance on searches for 
high-z spheroids 

Clearly, cosmic variance affects results based on small ar- 
eas such as the HDFs (see Sect. 4.1 and Figure 3). Daddi 
et al. (2000; hereafter DA) showed that the surface density 
variations of the Extremely Red Objects (EROs) popula- 
tion found in their survey can contribute to explain the dis- 
crepant results obtained by different authors about the den- 
sity of z > 1 ellipticals. Similarly, Eisenhardt et al. (2000), 
using a selection of z > 1 elliptical galaxies based on the 
(J — K) colour, examined the surface density variation be- 
tween their 7f -selected survey (124 arcmin 2 ) and the HDF- 
N and argued that clustering can affect results based on 
small area surveys. Actually, the cosmic variance is clearly 
apparent in our HDF-S and HDF-N catalogues (see Figure 
3, see also Lee & Hwang 2000). 

Some insights about the effects of cosmic variance on 
our conclusions can be obtained by comparing results based 
on our small (11 arcmin 2 ) area with data on the surface 
density of EROs derived by DA for a much larger sky area 
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Figure 7. Distributions of the density of stellar mass per unit comoving volume generated in the redshift bin as a function of redshift 
(in Mq Mpc~ 3 born in the redshift bin). The dashed lines represent the cumulative distributions of the mass density generated as a 
function of redshift, on a linear scale ranging from 0% to 100%. The dotted horizontal lines mark the 33 and 66 percentiles. 



(~ 700 arcmin 2 ). DA define as EROs those objects with 
colours (R — K) > 5. They could be passively evolving 
galaxies at z > 1 (with red colours because of the large 
K-correction) or heavily dust-reddened, star-forming galax- 
ies or AGN. Using the above definition of EROs, DA found 
a surface density of 0.67 EROs arcmin~ 2 at the limit of 
their survey (K < 19.2). A more stringent colour criterion 
of (R — K) > 6 should select objects at z > 1.3. At the 
limit of (R — K) = 6 they found a surface density of only 
0.10 EROs arcminT 2 . If we adopt the same magnitude limit 
of K < 19.2, the corresponding surface densities for z > 1 
and z > 1.3 ellipticals that we find in our sample are: 0.88 
arcmin" 2 (R — K > 5) and 0.17 arcmin~ 2 (R — K > 6). 
Thus, in spite of the very different survey areas and selec- 
tion criteria, the EROs statistics and our results about the 
fraction of E/S0 galaxies at high-z turn out to be in quite 
good agreement. This supports our conclusions about the 
demise of bright massive spheroids at high redshifts. 

This fast decrease of the areal density of E/S0 at z > 
1.5 is likely to correspond to a fast evolution of the mass 
function at these epochs. Altogether, a very active phase 
not only for the assembly of massive E/S0 galaxies in the 
field but also for the origin of stellar population seems to 
happen around z ~ 1 to 2. 



5 CONCLUSIONS 

We have analyzed a K-band flux-limited sample of 69 E/S0 
galaxies in three different areas (HDF-North, HDF-South 
and the NICMOS field, for a total of ~ 11 arcmin 2 ), where 
deep HST data are available. This sample is uniquely suited, 
as for the quality of the imaging and photometry, to the 
study of high-redshift spheroidal galaxies. 

Our main conclusions can be summarized as follows. 



• Wc confirm the finding by FA98 that massive E/SOs 
tend to disappear from flux limited samples at z > 1.4. At 
least the assembly of massive spheroids has to take place 
efficiently in the redshift interval 1 < z < 2. 

• The large spreads in the rest-frame colour distributions 
(both V — K and B — J) for the early-type galaxies observed 
at z ~ 1 in our sample indicate a wide range of ages, typi- 
cally from 1 to 5 Gyr. These objects are inconsistent with a 
very high redshift of formation for the bulk of stars, and are 
consistent with protracted star-formation (either continuous 
or episodic) down to z < 1. 

• Averaged over our (small) selected area, the bulk of 
stars (40%) are formed in the redshift interval 1 < z < 2. Of 
the remaining fraction, half form at z < 1 and half at z > 2. 

• Our finding of a demise of E/SOs at z > 1.5 in small 
fields is confirmed on much larger areas by a similarly fast 
decrease in the number density of EROs, when the colour 
limits are changed from (R — K) > 5 to (R — K) > 6 (cor- 
responding to an increase of the redshift cutoff from z ~ 1 
to z ~ 1.3, for typical evolved stellar spectra). 

• Hierarchical models for the formation of spheroidal 
galaxies may naturally explain the present results. 
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Table 1. Data on the sample galaxies 
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